arXiv:1508.06740vl [astro-ph.CO] 27 Aug 2015 


Primordial non-Gaussianities after 
Planck 2015: an introductory review 

Sebastien Renaux-Petel 1,2 

1 Institut d’Astrophysique de Paris, UMR-7095 du CNRS, Universite Pierre et Marie Curie, 
98 bis bd Arago, 75014 Paris, France. 

2 Sorbonne Universites, Institut Lagrange de Paris, 98 bis bd Arago, 75014 Paris, France 

Abstract. Deviations from Gaussian statistics of the cosmological density fluctuations, so-called 
primordial non-Gaussianities (NG), are one of the most informative fingerprints of the origin 
of structures in the universe. Indeed, they can probe physics at energy scales inaccessible to 
laboratory experiments, and are sensitive to the interactions of the field(s) that generated the 
primordial fluctuations, contrary to the Gaussian linear theory. As a result, they can discriminate 
between inflationary models that are otherwise almost indistinguishable. In this short review, 
we explain how to compute the non-Gaussian properties in any inflationary scenario. We review 
the theoretical predictions of several important classes of models. We then describe the ways NG 
can be probed observationally, and we highlight the recent constraints from the Planck mission, 
as well as their implications. We finally identify well motivated theoretical targets for future 
experiments and discuss observational prospects. 

Les deviations a la gaussianite des fluctuations cosmologiques de densite, ou non-gaussianites 
primordiales, nous fournissent des indices precieux quant a l’origine des grandes structures de 
l’univers. Elies permettent en effet de sonder la physique a des echelles d’energie inaccessibles 
en laboratoire, et sont sensibles aux interactions du (ou des) champ(s) a l’origine des fluctua¬ 
tions primordiales, contrairement a la theorie lineaire gaussienne. Elies nous permettent ainsi de 
differencier des modeles autrement presque indistinguables. Dans cette courte revue, nous ex- 
pliquons comment calculer les proprietes non-gaussiennes des fluctuations generees pendant tout 
scenario d’inflation. Nous passons en revue les differentes predictions theoriques de plusieurs 
grandes classes de modeles. Nous decrivons ensuite la fagon dont les non-gaussianites peuvent 
etre contraintes observationnellement et nous soulignons les contraintes recentes apportees par 
la mission Planck , ainsi que leurs implications. Nous discutons enfin les perspectives observa- 
tionnelles en identifiant des objectifs realistes et motives theoriquement. 
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1 Introduction 

Thanks to unprecedented observational efforts in the last two decades, we now have at hand high 
quality data of the two main cosmological probes, namely the Cosmic Microwave Background 
(CMB) anisotropies and the Large Scale Structures (LSS). The picture of the primordial universe 
that emerges from these data is surprisingly simple: all current observations are consistent with 
the A-Cold Dark Matter model, with initial conditions provided by the simplest inflationary 
models. The only ingredient of the latter is a single canonical scalar field minimally coupled 
to gravity and evolving on top of a very flat potential (we call them single-field slow-roll in the 
following). These scenarios provide a very good fit to the data (see the contribution by Martin, 
Ringeval and Vennin to this volume) while alternatives to the inflationary paradigm are less 
compelling (see the contribution by Peter and Lilley). However, despite being phenomenologi¬ 
cally successful, these models can not ultimately be considered as satisfactory, for at least two 
reasons: they are decoupled from the rest of physics, and they lack a ultraviolet completion. 

In particular, as soon as one wants to embed the inflationary paradigm into quantum field 
theory, it becomes surprisingly challenging to realise slow-roll models in an honest-to-God way 
[1] (see also the contribution by Silverstein to this volume). One is then led to consider more 
complicated inflationary models, involving for instance several degrees of freedom, non-canonical 
actions or features in the potential. While these models can easily be made degenerate with 
the simplest ones at the leading-order approximation, their complexities often reveal themselves 
at next-to-leading-order in the form of so-called primordial non-Gaussianities: deviations from 
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Gaussian statistics of the primordial density fluctuations. Cosmological data now put stringent 
bounds on them, constraining inflationary models in a way that would be otherwise impossible. 
Here we explain how to compute the non-Gaussian properties of inflationary models and we 
review the associated theoretical predictions in several important classes of models. We mention 
the ways non-Gaussianities can be probed observationaly, as well as the current and expected 
future constraints. This field of research has been very active in the past decade and we can not 
pay entire justice to it in this short introductory review. For complementary details, we refer the 
interested reader to the reviews [2-12], and to the references to the original literature therein. 
We use units in which c = h = Mp\ = 1. 

2 Inflation and the origin of the large scale structure of the universe 

From quantum fluctuations to primordial perturbations. Before discussing primordial 
non-Gaussianites (NG), it is appropriate to remind the reader of some basic facts about the 
inflationary origin of the large scale structure of the universe, referring her/him to the other 
contributions in this volume for a more detailed treatment. On cosmological scales, the geometry 
of the universe can be described at first approximation by a homogeneous and isotropic metric 
of the (flat) Friedmann-Lemaitre-Robertson-Walker (FLRW) form: 

ds 2 = —d t 2 + a{t) 2 8ijdx l dx 3 (2.1) 

where t is the cosmic time, a denotes the scale factor of the universe, and x l are spatial comoving 
coordinates. The large scale structure of the universe is then described by small fluctuations 
above this background, which we label by their fixed 3d Fourier wavevectors k in comoving 
space. The conceptual problems of the hot Big-Bang model all derive from the increase of the 
comoving Hubble radius ( ciH ) _1 , where H = a/a denotes the Hubble scale. Inflation solves them 
by providing an earlier phase in which ( aH ) _1 decreases sufficiently that the observable universe 
was initially inside the Hubble radius (see Fig. 1 and the contribution by Ellis and Uzan to this 
volume about the causal structure of inflation). This way, scales of cosmological interest have 
their physical wavelength A = less than H at the beginning of inflation, exit the Hubble 
radius during inflation before re-enteringjt during the Big-Bang era. In this process, unavoid¬ 
able sub-Hubble quantum fluctuations <5</(. of the field(s) driving inflation become classical after 
Hubble exit [13], manifesting themselves in the Big-Bang eras as stochastic fluctuations of the 
number density nx of the various components of the universe: radiation, neutrinos, cold dark 
matter, baryons ... Remarkably, all curent data indicate that there is no spatial fluctuation of 
the composition of the primordial plasma on cosmological scales, i.e. 5 = 0 for X,Y. 

This fact, referred to as the adiabaticity of the primordial fluctuations, implies that the latter 
can be described by one only quantity: the so-called curvature perturbation on uniform density 
hypersurfaces ( deep in the radiation era. It also points to a common origin of these fluctua¬ 
tions, like in the simplest models of inflation driven by a single scalar field. In these scenarios, ( 
becomes constant soon after Hubble exit during inflation, which enables one to make predictions 
at a much later time (and lower energy) in the radiation era without having to deal with the 
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Figure 1. The inflationary expansion turns sub-Hubble quantum fluctuations into classical super-Hubble 
perturbations, whose effects we observe in the CMB and the LSS. 


complicated dynamics of the intermediate phase of reheating after inflation. As we will see, the 
super-Hubble conservation of the curvature perturbation is lost in general when more than one 
degree of freedom is active during inflation, like in multifield inflationary models. One should 
then keep track of cosmological fluctuations on super-Hubble scales until an adiabatic limit is 
reached, which requires more attention. 

Background dynamics of inflation. — An expanding universe with shrinking comoving Hub¬ 
ble radius is accelerating: a > 0. In such a phase, the dimensionless parameter e = —H/H 2 
is less than 1. For simplicity, and in agreement with cosmological data, one assumes in the 
following that the inflationary expansion is quasi-exponential, i.e. that the Hubble parameter 
varies sufficiently slowly that e< 1 and r/ = e/(He) <C 1. More generally, we assume except 
otherwise stated that all quantities during inflation evolve much less rapidly than the scale factor: 
X/X < a/a = H. At leading-order in this slow-varying approximation, H ~ cst and a(t) oc e Ht . 
In terms of the conformal time r defined by d t = adT, this gives a(r) ~ —1 /(Hr), where r 
grows from —oo to 0 during inflation. Additionally, in terms of the so-called number of e-folds N 
defined such that diV = H df, or equivalently a oc e N , the resolution of the conceptual problems 
of the Big-Bang model requires that N > 60. At the level of classical field theory at least, it 
is straightforward to achieve a phase with the aforementioned properties by using a scalar field 
rolling down a sufficiently flat potential. However, there are many other proposals to achieve 
such a phase and in the following, we show how studying perturbations on top of this background 
helps to differentiate amongst competing models. 
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3 Primordial power spectrum and higher-order correlators 

Cosmological observations enable one to constrain the large-scale statistical properties of the 
curvature perturbation. Its two-point correlation function in Fourier space defines the primordial 
power spectrum P((k) such that 

(Cki Ck 2 ) = (2vr) 3 i5(k 1 + k 2 ) P((k) (3.1) 

where k = |ki| = |k 2 1- On cosmological scales 0.008 Mpc^ 1 < k < 0.1 Mpc -1 , it can be accurately 
described by a power law 

- 2 P d k) - Mk.) (-) (3.2) 

with the pivot scale &:* = 0.05 Mpc^ 1 and the constraint [14] 

n s = 0.968 ± 0.006 (68%CL). (3.3) 

This convincing detection of a percent-level deviation from perfect scale invariance of the primor¬ 
dial power spectrum (n s = 1) is one of the most important recent achievement in observational 
cosmology. It agrees with the predictions of the simplest model of inflation and, together with 
constraints on the amplitude of primordial gravitational waves, it enables one to pin down the 
range of inflationary scenarios responsible for the generation of these fluctuations (see the con¬ 
tribution by Martin, Ringeval and Vennin to this volume). 

In addition to this, one can extract complementary information on the interactions of 
the field(s) driving inflation, through the use of higher-order connected correlation functions of 
the curvature perturbation — which would vanish for perfect Gaussian statistics — like the 
bispectrum 

(Ckj Ck 2 Ck 3 > = (2 vt) 3 <5(^ k;) B c ( ki) (3.4) 

i 

and the trispectrum 

(Cki Ck 2 Ck 3 Ck 4 )c = (27 t) 3 «5(^ kj) T c (ki). (3.5) 

i 

The factors <5(^Tki) above result from statistical homogeneity and imply that the wavevectors 
k i respectively sustain a triangle and a tetrahedron in Fourier space. Additionally, statistical 
isotropy entails that the orientations of these polyhedra are irrelevant, so that only their shapes 
and overall scales matter. This still leaves 3 independent variables to describe the bispectrum 
(and 6 for the trispectrum), which highlights the richness of information contained in these cor¬ 
relation functions compared to the power spectrum (3.1), that depends on only one variable. 

Now concentrating on the bispectrum, it is customary to define its profile S by 

(3.6) 
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where A s denotes the amplitude of the curvature power spectrum (see Eq. (3.2)). The overall 
magnitude of the dimensionless function S, which is denoted as the non-linearity parameter 
/nLj provides an estimate of the importance of the bispectrum. The simplest models of inflation 
generate a tiny amount of non-Gaussianities [15-17], with /nl = 0(e, rj) <C 1, under the following 
hypotheses: a single field drives inflation and generates the primordial fluctuations; it has a 
standard kinetic term and a smooth potential; its perturbations are in a vacuum state. However, 
as we will see, violating any one of these restrictive conditions can result in observably large NG, 
|/nl| 1, while recent constraints give the interesting bound /nl 0(10). 

Besides this qualitative estimate, the detailed geometrical dependence of S on the con¬ 
figuration of the wavevectors enables one to differentiate amongst the different classes of in¬ 
flationary models beyond single-field slow-roll [18]. For example, scenarios with multiple light 
degrees of freedom generate a bispectrum than can be large in squeezed configurations, i.e. for 
Aq Aq — Aq, in contrast with single-field models, which predict a vanishingly small bispectrum 
in this limit. Additionally, models with important derivative interactions, for instance with non- 
canonical kinetic terms, are characterized by a bispectrum that is maximum for wavevectors of 
similar magnitude, i.e. around equilateral configurations Aq ~ Aq — Aq. Scenarios with features 
also come hand with hand with their signatures in the form of specific oscillatory bispectra. This 
way, cosmologists have identified in the past decade a useful dictionary between several broad 
classes of inflationary mechanisms for the generation of primordial fluctuations and the associated 
type of non-Gaussianities. In this respect, it is useful to decompose the geometrical dependence 
of S in two parts: its shape and its scale-dependence. In the simplest ‘smooth’ models of infla¬ 
tion, without features, modulations of the potential or excited states for instance, S/Aq, Aq, Aq) 
depends mainly on the ratios between the norms of the wavectors, e.g. Aq/Aq and Aq/Aq, and 
mildly on their overall scales. In the following, we order the A:/s such that Aq < Aq < Aq, so that 
this shape information can be represented by plotting the two-dimensional function 5(1, X 2 , £ 3 ), 
where 0 < £ 2,£3 < 1 and 1 < £2 + £3 (to satisfy the triangle inequality). In more complicated, 
non-scale invariant, models, this has to be complemented by the study of the dependence of S 
on the overall scale K = Aq + Aq + Aq while keeping the ratios Aq/Aq and Aq/Aq fixed. 

4 Methods to calculate primordial non-Gaussianities 

In this section, we explain how to calculate higher-order correlation functions in any inflationary 
scenario. For definiteness, we consider models involving scalar fields and a cosmological metric 
that takes the flat FLRW form at the background level. However, the methods we present are 
generic and are easily applicable to models with other matter contents, such as vector fields, and 
other background symmetries. 

Cosmological perturbation theory. — We start by decomposing the metric and the fields < -p 1 
(I = 1, ...,M) into their background and fluctuating parts: 

g^(,t,x l ) = g^{t) + 5g^{t,x l ) 

</> 7 (A, £*) = </ 7 (t) + <h/ 7 (t, £*) 
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(4.1) 

(4.2) 


where the background parts only depend on cosmic time. We then expand the (gravitational 
and matter) action about the background solutions of the equations of motion, which we write 
schematically as: 


S 


S + sWtfg^M 1 ) + S (3) (%„ H 1 ) + S^(Sg^, 5</>*) + .. 

'-V- 


(4.3) 


The second-order part governs the behaviour of linearised fluctuations, while the higher-order 
parts encode their interactions. A crucial subtlety in these calculations comes from the ambiguity 
in identifying spacetime points between the idealised background and the true perturbed uni¬ 
verse, and the gauge freedom it entails [19, 20]. As a consequence, one must pay special attention 
to the proper identification of the true dynamical degrees of freedom and the construction of the 
corresponding gauge-invariant variables. For instance, a model with canonical Einstein-Hilbert 
action and N scalar fields contains 4 + N (Lagrangian) degrees of freedom: two tensor modes, 
so-called gravitational waves, two tensor modes, and N scalar modes. These different sectors are 
decoupled at the linearised level and evolve very differently. Tensor modes are only sensitive to 
the background spacetime geometry, i.e. to the evolution of the scale factor. Without further 
sources, vector modes decay due to the expansion of the universe. On the contrary, details of 
the evolution of the scalar perturbations depend on the model at hand. They provide the main 
sources of the curvature perturbation (, which imprinted the CMB anisotropies and seeded the 
formation of LSS. We concentrate on them in the following, and discard the other types of fluc¬ 
tuations for simplicity of presentation. 


Quantisation. - Once the second-order action is known, one can proceed to its canonical quan¬ 
tisation by promoting the propagating degrees of freedom to quantum operators. The procedure 
is well known and is straightforward in the single-field case [21] but it should be taken care of 
carefully in multifield scenarios (see e.g. Ref. [22]). The normalisation of fluctuations, which 
is arbitrary in a classical context, is fixed upon imposing the canonical commutation relations. 
Eventually, a choice of vacuum |0) should also be made: deep inside the Hubble radius (k 3> aH), 
fluctuations are not sensitive to the expansion of the universe and, following the equivalence 
principle, it is legitimate to favor the approximate Minkowski vacuum, also called Bunch-Davies 
vacuum [23] (see below for a discussion of the impact of an excited initial state). This completes 
the quantisation procedure, leaving only the task of accordingly solving the linear equations of 
motion deduced from the second-order action, either analytically or numerically. Finally iden¬ 
tifying the vacuum expectation value of observables in the quantum theory with the statistical 
ensemble average of the corresponding variables in the classical theory, one can then extract 
predictions, for the curvature power spectrum (3.1) in particular. 

In-in formalism. After having quantified the linear Gaussian theory, and identified the 
interacting action in Eq. (4.3), one can determine higher-order correlation functions using the 
so-called in-in (also called Schwinger-Keldysh) formalism [24-26]. Starting from first principles in 
quantum field theory, it shows that the expectation value of an observable 0(t) can be computed 
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perturbatively as 


(in\0(t)\in) = 


T exp i 


' — oo(l— ie) 


H^dt' 


o ! {t) 


T exp —i 


' —oo(l -\-ie) 


H T {t")dt h 


| 0 ) 


where | in) is the vacuum of the interacting theory at some moment ti in the far past, T denotes 
the time-ordered product, the J’s indicate the use of the interaction picture and Hj is the 
interacting Hamiltonian. At first-order in the latter, as relevant for the calculation of the tree- 
level bispectrum, one finds 


(0(t)} (1) = 2 Re 


-i / df / (0|O / (£).£//(f , )|0) 


whereas this reads 


(0(t))W = d t' dt"(0\H I (t')O I (t)H I (t")\0) 


' —oo J —oo 


- 2 Re 


df' / dt ,, (o|o / (t)fr J (t , )fr / (t ,, )|o) 


' —oo J — oo 


(4.4) 


at second order, which is necessary for computing the trispectrum for instance (see Ref. [27] for 
general expressions at higher orders). These terms are then simply evaluated by applying Wick’s 
theorem (let us recall that all fields here are in the interaction scheme and hence are free fields). 
Eventually, note that, as usual in quantum field theory [28, 29], as one wishes to compute these 
expectation values in the vacuum of the interacting theory, and not in the vacuum of the free 
(linear) theory, one should slightly deform the integration contours in the complex plane (the 
standard ie prescription), which turns off the interactions in the far past and renders the time 
integrals well-behaved then. 


5N formalism.— The in-in formalism is very general and can be a priori applied to any 
inflationary model. However, in multifield models, it is often difficult to obtain accurate analytical 
expressions for the linear mode-functions beyond a few e-folds after the time of Hubble crossing. 
For this reason, it is useful to have at hand another formalism called 5N [30-34], The latter 
relies on the separate-universe picture [35]. It states that each super-Hubble patch — as relevant 
for dealing with fluctuations on super Hubble-scales k <C aH — evolves like a separate FLRW 
universe which is locally homogeneous and evolves independently from its neighbours. However, 
patching these regions together enables one to track the evolution of the curvature perturbation 
on large scales, just by using background quantities. More specifically, one can express the 
curvature perturbation £(t, a;*) as the difference between the number of e-folds of expansion N 
between an initial flat arbitrary hypersurface and the uniform energy density hypersurface at 
time t, and its corresponding background value: 

£(£, x l ) = N(t , £*; x l ) - N(t , £*). (4.5) 
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This formula is completely generic and non-perturbative. Let us now specify it to attractor 
solutions of multifield inflationary models 1 . Their dynamics is entirely dictated by specifying 
the initial values of the scalar fields, i.e. the spatial dependence in N(t,t*;x‘ l ) can be encoded 
as N(t, One can then Taylor-expand Eq. (4.5) as a function of the field fluctuations 

on the initial flat hypersurface. By writing (f> A = </>* + Q A , one obtains 

c = N A Q A + \n ab Q a Q b + \n abc Q a Q b Q c + ... (4.6) 

2 6 


where 


N a = 


d N 


d(p A 


N ab 


d 2 N 
d(p+ dd)+ 



(4.7) 


From Eq. (4.6) it is straightforward to deduce the statistical properties of the curvature pertur¬ 
bation in terms of the ones of the scalar fields. Defining their spectra, bispectra and trispectra 2 : 


(QkQvY = C AB (k ) (4.8) 

(Qi Ql QgY = B ABC {ki) (4.9) 

{Qt Qk 2 QZqO'c = T ABCD (k t ) (4.10) 


one obtains the expressions of the tree-level primordial spectrum, bispectrum and trispectrum 
as [37, 38] 

P ( (k) = N A N B C AB (k ) (4.11) 


B c (k 1 ,k 2 ,k 3 ) = N A N B N c B ABC (k 1 ,k 2 ,k 3 ) + N A N BC N D [C AC (k 2 ) + 2 perms.] (4.12) 
and 


T c (k*) = NaNbNcNdT^ 00 ^) + N AB N C N D N E [C ac ( h)B BDE (k 12 , k 3 , k 4 )\ 
+n ab n cd n e n f \C BD (k 13 )C AE (k 3 )C GF (k 4 )\ +n abc n d n e n f [C AD {k 2 )C BE {k 3 )C CF {k 4 )] 

-|-permutations (4-13) 

where kjj = |k ? ; + k ? | and ‘permutations’ indicate that one should add the respectively 11, 11 and 
3 inequivalent permutations of the kj to the last three lines, so that the trispectrum is totally 
symmetric. 

Interpretation.- - For definiteness, let us now specify the time of the initial flat hypersurface to 
be a few efolds after the period of Hubble crossing for the relevant cosmological scales. The in-in 
formalism then enables one to accurately predict the correlation functions of the field pertur¬ 
bations (4.8)-(4.10). The expressions above, although formal, yet deliver a useful information. 
For instance, one can see that the bispectrum (4.12) is a sum of two contributions: the first 
line, stemming from the linear term N A Q A in Eq. (4.6), results from the linear transfer to the 

1 A generalised SN formalism has been developed to deal with more complicated background dynamics [36]. 

2 A prime on correlation functions indicate that one omits the ever-present factor (27r) 3 (5(^ i ki). 







curvature bispectrum of the intrinsic bispectra of the field fluctuations B abc around Hubble 
crossing. The second line, proportional to the second-order coefficient Nab hr the expansion 
(4.6), is present even for perfectly Gaussian field fluctuations, and comes from the non-linear 
relation between the latter and the curvature perturbation £. A similar discussion holds for the 
trispectrum (4.13): the first line simply results from the linear transfer of the field trispectra to 
the curvature perturbation. The last two lines, on the opposite, solely come from the respectively 
quadratic and cubic terms in the expansion (4.6). The second contribution, of a type that is not 
present for the bispectrum, is more complicated, as it mixes the intrinsic bispectra of the field 
fluctuations and the nonlinear (quadratic) term in the relation (4.6). 

Barring it, one can therefore classify the contributions to the bispectrum and trispectrum 
in two main types. The first ones originate from purely quantum effects around (and possibly 
before) the time of Hubble crossing. Their precise geometrical forms depend on the model at 
hand, and their determinations require the use of the in-in formalism. They can be important 
in either single or multi-field scenarios, in models with large derivative interactions in particular 
[39-42], The second type comes from the classical non-linear relation between the curvature 
perturbation and the field fluctuations around the time of Hubble crossing. Their shapes are 
entirely fixed by the two-point correlation functions of the field fluctuations C AB , and therefore 
take universal forms in the simplest models, as we will see. Their importance relies on large 
non-linear terms in the expansion (4.6), and therefore a large sensitivity of the number of e-folds 
of expansion to the initial field values. No such sensitivity is present in (attractor) single-field 
models, so that this type of contribution is characteristic of multi-field models with non-trivial 
dynamics on super-Hubble scales (see below for details). 

5 General single-field models 

In this section, we discuss the non-Gaussian properties of the curvature perturbation generated in 
single field models of inflation. For simplicity, we restrict ourselves to the representative class of 
models known as ‘/c-inflation’ [43, 44]. They generalise the inflaton Lagrangian from its canonical 
form C = X — V((/)), where X = — ^g^d^dycj) is the field’s kinetic term, to a general function 
C = P{X,(f>). Such models can support a phase of inflation of a very different type from slow- 
roll scenarios. For example, the prototypical model of string-inspired Dirac-Born-Infeld (DBI) 
inflation [39, 45], with P = —1 /f{(j>) ^-y/l — 2 f(cj>)X — 1^ — V{4>) and a positive can achieve 
inflation on top of very steep potentials, not through Hubble friction like in standard models, 
but because the noncanonical kinetic term bounds the speed of the inflaton, irrespective of the 
potential: 2X = (j) 2 < l//(0). 

Although it may not seem intuitive at first sight, it is useful in single-field models to choose 
a coordinate system in which the field </> has no fluctuation 5<f>. It is always possible to choose 
such a coordinate system, called the uniform inflaton gauge, and to additionally impose that the 
spatial part of the metric takes the form (still discarding tensor perturbations) 

9ij = a 2 (t)e 2C(t,3:l) 5ij . (5.1) 
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The variable £ defined in this way coincides in this gauge with the gauge-invariant curvature 
perturbation we are interested in. It is therefore quite convenient to directly deal with it in 
calculations, without having to make explicit its relation with 5<j> at non-linear orders. By 
solving the constraint equations to express the non-dynamical parts of the metric in terms of the 
propagating degree of freedom £ (this is most easily done by using the ADM formalism [46]), one 
can derive the second-order action 



(5.2) 


where 

2 = p ,x 

p x + 2 XP )X x 


(5.3) 


is called the (squared) sound speed of perturbations. It equals one in canonical models, in 
which P = X — V((pi), but it can be much less than one in general, like in the example of DBI 


models where c s = y 1 — f(cj))(j) 2 . At leading-order in the slow-varying approximation, one can 
analytically solve the equation of motion deduced from Eq. (5.2): 


Cfe(r) = 


iH * 


\/4e*c s */c 3 


(1 + ikc s *r) e 


—ikc s *T 


(5.4) 


where * denotes evaluation at the time of the sound horizon crossing kc s = aH. As explained 
above, one has chosen Bunch-Davies initial conditions and the normalisation comes from the 
quantization procedure. It is then straightforward to deduce the primordial power spectrum 

k 3 k 3 / u2 




(5.5) 


with a mild scale dependence given by the slight dependence of H, e and c s on the time of eval¬ 
uation *. 


Equilateral and orthogonal non-Gaussianities. — Now moving on the bispectrum, one first 
calculates the leading-order cubic action [40, 47]: 


c( 3 ) _ 

d l.o. — 


dtd^x a A e ( -^ — 1 


- V 

a c s 


+ / dt d 3 x cT 


He 2 


where one concentrated on the main sources of non-Gaussianities, 

£ = XP x + 2 X 2 P xx = H ^ L 


\ = X 2 Pxx + -X s P X xx 


^ - 1 - Y ) C ' ’ (5 - 6) 


(5.7) 

(5.8) 


and all parameters are taken to be constant in this approximation. Using Eq. (4.4), one can then 
deduce the shape of the bispectrum as S = S x + S c , with 


S A = ^ - 1 - — 


2A\ 3k±k2ks 


2 K 3 


(5.9) 
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and 


S c = -s - 1 


1 


k\k 2 k 3 


JC £ + 2I^Y, k itf + l E k i 

i>j v£j i 


(5.10) 


and where K = k\+k 2 +k 3 . One therefore hnds /nl = O ^ — 1, , i.e. large non-Gaussianities 



(a) Shapes of S x in Eq. (5.9) (left) and S c in Eq. (5.10) (right). 



(b) Shapes of S eq in Eq. (5.11) (left) and of the absolute value of 5 orth in Eq. (5.12) (right). 

Figure 2. Different shapes S x (1, x 2 , x 3 ) as a function of (x 2 , x 3 ). We set them to zero outside the region 
1 — x 2 < x 3 < x 2 and normalize them to one for equilateral triangles x 2 = x 3 = 1. 


/nl 1 are generated when the deviations to the non-canonical structure P = X — V(4>) are 
significant. The two shapes S x and S c are represented on the top panel of Fig. 2, where one can 
see that they are similar, both peaking in the equilateral limit. In a first approximation, one can 
represent both of them by a unique simpler template called equilateral [48] (and whose form is 
motivated by data analysis, see section 8): 
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Although this is not obvious from the cubic action (5.6), one can use the linear equation of motion 
deduced from the second-order action to express in terms of the interactions £ 3 and C{diC) 2 
only [49] (this becomes more transparent in a gauge where the scalar field itself is perturbed, 
and one can neglect the mixing with gravity for large NG). These derivative interactions are 
suppressed when any individual mode is far outside the Hubble radius, either by gradient terms, 
or because of the super-Hubble conservation of £). This should not come as a surprise then that 
the bispectrum is maximal for modes that cross the sound horizon at approximately the same 
time, i.e. for equilateral triangles. 

By looking at Fig. 2, one can see that there are small differences between the shapes S x and 
S c , especially around flattened triangles fa + fa ~ A; i. One can therefore highlight their differences 
by considering an appropriate linear combination of them that subtracts their similarities. This 
way, one can cover more efficiently than with S eq alone the two-dimensional space of shapes 
spanned by S x and S c . A simple template for such a shape, called orthogonal, is given by [50]: 


27 

5 orth 1 .forth 

-|^q J NL 


-( 


*? 


\fafa 


+ 2 perm. + 


fa \ 8 ' 

--h 5 perm. - - 

k 2 J 3 


(5.12) 


The two shapes S eq and S ,orth are represented on the bottom panel of Fig. 2 (we actually plot 
the absolute value of S' orth so that its difference with the equilateral ansatz for flattened trian¬ 
gles is more visible). Although we concentrated on ^-inflationary scenarios of the type P(X,(j>), 
it should be clear from the above discussion that these shapes are signatures of derivative in¬ 
teractions in general. They indeed emerge in more general higher-derivative scenarios, such as 
ghost-inflation [51], Galileon-like models [52-54], or Horndeski and generalised theories [55-57]. 


Consistency relations. — We have focused above on the characteristic equilateral-type shapes 
that single-field inflation with derivative interactions can generate. It is also interesting to con¬ 
sider what type of non-Gaussian signal single-field inflation cannot generate. In this respect, one 
can show that, under very general conditions, a phase of inflation driven by a single scalar field 
generates a bispectrum such that [58] 


lim fa, fa) = (1 - n s (ki))P ( ;(k 1 )P(;(fa) (5.13) 

Ay —>0 


where 


n s (k) - 1 


din [k 3 P^(k)] 
d In k 


(5.14) 


is the scalar spectral index. In other words, the squeezed limit of the bispectrum is suppressed by 
1 — n s and vanishes for a perfectly scale-invariant power-spectrum. Given the observed smallness 
of the deviation from scale-invariance (see Eq. (3.3)), this implies that a detection of a large 
bispectrum signal in the squeezed limit (f^\ j > 1) would rule out all models of inflation based 
on a single scalar field, irrespective of any details like the form of the potential or the kinetic 
term! Given the importance of this relation, let us give a sketch of its proof. It relies on the 
fact that a super-Hubble conserved curvature perturbation locally acts as a background field, 
only rescaling the spatial coordinates within a given Hubble patch (see Eq. (5.1)). By writing 
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(CkiCk 2 Ck 3 ) = ((CkiCk 2 )Ck 3 Ck 3 ) ) where (CkiCk 2 )Ck 3 is the expectation value of CkiCk 2 given that 
Ck 3 has a particular value, one can see that the squeezed limit bispectrum (k-^ —> 0) measures 
the effect of a long-wavelength fluctuation on the short wavelength power spectrum. From the 
argument above, this effect is proportional to the change of the power-spectrum under a rescaling 
of the spatial coordinates, or equivalently, to its deviation from perfect scale-invariance. 

As should be clear from this sketch, its crucial argument is the constancy of £ on super- 
Hubble scales. In this respect, it should be said that, strictly speaking, single-field inflation only 
implies the existence of a super-Hubble constant mode for its other mode usually decaying 
exponentially. However, by choosing ad hoc potentials and initial conditions, it is possible to 
render this decaying mode significant on super-Hubble scales, thus violating the consistency 
relation (see e.g. [59, 60]). This mild limitation does not however limit the importance of the 
relation (5.13), which has been derived using various symmetry arguments, as well as generalised 
to higher-order correlation functions and slightly different situations (see e.g. Refs. [61-68]). 


6 Multi-field models 


In this section, we consider models characterised by the presence of several light scalar degrees 
of freedom during inflation. A prototypical exemple is the one of a collection of canonical 
scalar fields (p 1 minimally coupled to gravity and interacting through a generic potential V^cp 1 ). 
Contrary to the single field case, one can not remove here all the fluctuations of the scalar fields 
by a suitable choice of coordinates. Instead, it is convenient in that case to use the so-called 
spatially flat gauge, in which the spatial part of the metric g tJ takes its background form and the 
field fluctuations Q 1 are made explicit. The second-order action then takes the form [30, 69-71] 


5(2) 


/ 


df d 3 x a 3 



(diQ z ) 2 \ 

a 2 ) 


f 1 

[a 3 ■ . I 


~^9l9J 


M u 



( 6 . 1 ) 


When the effective mass matrix Mjj is negligible compared to H 2 around Hubble crossing, i.e. 
in the case of effectively light degrees of freedom, it follows from Eq. (6.1) that the fluctuations 
Q 1 acquire independent almost scale-invariant power spectra soon after Hubble crossing: 


C AB {k ) 


Hi 

2k 3< 


xAB 


( 6 . 2 ) 


At this stage, IV-field inflation resembles N copies of single-field inflation. The crucial difference 
lies in the super-Hubble regime. One can indeed show that, for any relativistic theory of gravity 
and any matter content, the local energy conservation implies that, for linear perturbations about 
an FLRW metric [35]: 


C = -H^^ + O 
P + P 



(6.3) 
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Figure 3. In multifield inflation, fluctuations are conveniently decomposed into an instantaneous adia¬ 
batic mode, along the background trajectory, and instantaneous entropic/isocurvature modes, orthogonal 
to it [72]. 

where p and p are the matter energy density and pressure respectively, and 6p nar \ is the non- 
adiabatic pressure perturbation, defined as 

<5y>nad = dp- ^rdp. (6.4) 

P 

In (attractor) single-field models, all variables x share the same time shift, 6t = along a 
single phase-space trajectory. As a consequence, the so-called entropy perturbation between any 
two quantities, T xy = ^ — ^, vanishes, and so does the non-adiabatic pressure perturbation on 
large scales in particular. On the contrary, in a multifield model, there is a family of inflationary 
trajectories, and perturbations in directions orthogonal to the classical path in field space lead to 
cosmological evolutions that are not simply translations of it, sourcing the non-adiabatic pressure 
perturbation, and hence the curvature perturbation Q on super-Hubble scales [73]. 

A simple example is provided by two-field inflation. One can show in this context that 
^Pnad is proportional on super-Hubble scales to T^i^ = This quantity is proportional 

to the projection of the perturbations in the direction orthogonal to the background trajectory: 
Q s = — sin (6)Q 1 +cos (0)Q 2 , where one has introduced the angle 9 such that tan(0) = 4> 2 /tfi 1 (see 
Fig. 3). This type of perturbation is called entropic, or isocurvature. On the contrary, one can 
show that the curvature perturbation is given on large scales by £ = ——j====^=Q a: where 

Q a = cos(0)Q 1 +sin(0)Q 2 , called the adiabatic perturbation, is the projection of the perturbations 
along the direction of the background trajectory. In details, the equation (6.3) reads, in this 
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context [72]: 


c = —2 


ho 


hl\2 


+ 


Qs + O 


2\2 



(6.5) 


In other words, the isocurvature perturbation sources the adiabatic curvature perturbation on 
super-Hubble scales as soon as the trajectory bends in field space, i.e. 0 ^ 0. This genuine mul¬ 
tifield effect can drastically modify the properties of the cosmological perturbations and should 
be accounted for in any model with multiple degrees of freedom. 


Local non-Gaussianities. — Our discussion here was at the linear level up to now but can be 
extended to non-linear orders (see e.g. Refs. [74-76]): NG of the curvature perturbation can be 
generated by the non-linearities of the transfer mechanism from the super-Hubble isocurvature 
perturbations, or can be linearly inherited from intrinsic non-Gaussianities in the isocurvature 
sector, which is less constrained than the adiabatic one by the slow-roll requirements. We used 
the bending of the trajectory in multifield inflation as a representative concrete example of this 
transfer mechanism [77-83], but it can arise at the end of inflation [84-87], or even after, like 
in the curvaton scenario [88-92]. However, whatever the specific model that is considered, the 
shape of the associated non-Gaussianities takes a universal form. Indeed, these processes arise on 
super-Hubble scales, where gradients can be neglected by definition, so that the physics become 
local in real space, and hence non-local in Fourier space, i.e. the bispectrum correlates large 
and small scale Fourier modes. This becomes transparent in the language of the 5N formalism 
presented in section 4, where these scenarios correspond to large super-Hubble classical non- 
linearities, i.e. to the second contribution to the bispectrum in Eq. (4.12). With the result (6.2) 
for the spectra of light scalar fields, this translates into the so-called local shape [93-95] (see 
Fig. 4) 

s ‘“ = + 2 perm ) m 


where 


, oc _ 5 N A N B N AB 
NL ~ 6 ( N c NC ) 2 


(6.7) 


Large non-Gaussianities can hence be generated in the squeezed limit -C ^2 — k±, where the 
shape (6.6) peaks, in sharp contrast with single-field models, in which the consistency relation 
(5.13) severely bounds the amplitude of the non-Gaussian signal in this limit. 

A similar discussion holds at higher orders in perturbation theory. In particular, the last 
two contributions in Eq. (4.13) generate in these models a trispectrum of the form 


T ; i oc (kj) = t nl [P c (fci 3 )P c (fc 3 )P c (A; 4 ) + perm.] + — g NL [Pt^P^k^P^k^ + perm.] (6.8) 


with 
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Figure 4. Shape of S loc in Eq. (6.6), with the same specifications as in Fig. 2. 


and 


9 NL = 


25 N A bcN a N b N c 


54 ( N d N d ) 3 ' ^ 6 ' 10 ' ) 

Similarly to the case of the local bispectrum, the tnl and ^nl shapes peak in particular ‘soft’ lim¬ 
its of the tetrahedron, with respectively a diagonal (internal wavevector) and one side of the gen¬ 
eral wavevectors being much smaller than the others. The parameter ^nl corresponds to a contact 
interaction and is related to intrinsic cubic non-linearities (the third term in Eq. (4.6)), while 
tnl corresponds to a scalar-exchange interaction and is related to second-order non-linearities, 
like /nl- The two are actually related by the inequality [96] 

2 


6 


tnl > ( ^/nl 


( 6 . 11 ) 


which is a particular case of a more general consistency relation [64]. The bound (6.11) is satu¬ 
rated in scenarios where a single field beyond the inflaton generates the curvature perturbation, 
and larger values indicate a more complicated dynamics. 


7 Beyond the simplest non-Gaussian shapes 

The non-Gaussian shapes we have described until now — local, equilateral, orthogonal — are 
well motivated and cover interesting broad classes of inflationary models. However, they only 
constitute the tip of the non-Gaussian iceberg. In this section, we briefly mention some other 
non-Gaussian signatures. 


Excited initial states. — In models with excited {i.e. non Bunch-Davies) initial states, the 
usual mode function (5.4) in e~ lkr acquires a small negative energy component ~ e lkr . This re¬ 
sults in (so far unobserved) oscillations in the primordial power spectrum [97-99] , which severely 
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constrains the magnitude of the effect. As for the bispectrum (4.4), its computation involves an 
integration of the product of three mode functions with wavenumbers k \, k 2 and k 3 . Hence the 
main correction to the Bunch-Davies result is to replace one of the ki s with —ki in this calcula¬ 
tion, turning the usual factors of 1/A' coming from the integration of e~ lKr into l/(k 2 + k 3 — k\) 
and its permutations. As a result, although detailed predictions are model-dependent (see e.g. 
[40, 100, 101]), the most important feature of this type of modification is to enhance the bis¬ 
pectrum for flattened triangles, with k 2 + k 3 ~ k\. A simple ansatz that captures this main 
characteristic is given by [ 101 ] 


gflclt 




+ 2 perm. 


(k 1 _ \ „ 

--1- 5 perm. + 3 

U’2 J 


(7.1) 


which is simply a linear combination of the equilateral and orthogonal templates. 


Scale-dependent resonance and feature models. — Periodic features in the inflationary 
potential naturally emerge when trying to UV-complete large field inflationary models, in par¬ 
ticular in a string theory context [102, 103]. These features can resonate with the sub-Hubble 
oscillations of the inflationary perturbations, generating an oscillatory bispectrum whose shape 
takes the form [104, 105] 

S' res = sin [ C hr (A* + k 2 + k 3 ) + 0 ] (7.2) 

in the simplest model, where C depends on the periodicity of the features. 

Temporary violations of the slow-varying evolution can occur if there are sharp features in 
the inflationary potential or sharp turns in field space in multifield inflation for example [106- 
109]. This can generate oscillatory correlation functions of a type different from the one in 
resonance models, the simplest bispectrum signal being of the form 

S icat = ~tit sin [w {ki + k 2 + k 3 ) + </>], (7.3) 

where oj again depends on the specific disruption of the attractor regime. Contrary to the simple 
almost scale-invariant shapes discussed in the previous sections, the two templates (7.2) and 
(7.3) provide examples of well motivated shapes whose scale-dependence is non-trivial. 


Quasi-single-field inflation. Technically natural models in supersymmetry often have extra 
fields beyond the inflaton with masses m of order the Hubble scale during inflation [27, 110— 
112]. In this sense, they interpolate between truly multifield inflation, where all fields are almost 
massless (m/H < 1), and single field inflation, which can be thought of as the limit m/H —> 00 . 
In these models, the conversion of non-linearities from the iscocurvature sector to the curvature 
one results in a bispectrum shape 


Sqsi 


3H/ fQS i (kj + kl + kl)(k 1 k 2 k^- l/ 

10 NL {h + k 2 + k 3 )t 3p 


(7.4) 
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that behaves in the squeezed limit as Sqsi oc , where v = y This scaling be¬ 

havior is intermediate between that of the local shape (k\/k .$) and that of the equilateral shape 
(k^/ki), exemplifying how soft limits can act as a particle detector. 

The trispectrum beyond the local shapes. — Beyond the local trispectrum (6.8) that arises 
in multifield inflation, the intrinsic quantum trispectrum (the first line in Eq. (4.13)) can be large 
in models with derivative interactions for instance [113-118], with shapes that are counterparts 
of the equilateral and orthogonal bispectra [119]. Additionally, classical and quantum effects 
combine in multifield models with derivative interactions (see the second line in Eq. (4.13)). 
Some of these models can hence generate a superposition of the local and equilateral bispectra, 
with a particular shape of the trispectrum whose amplitude is given by the product /^l/nl [120]. 

Inflation with gauge fields. A coupling of the inflaton to the kinetic term of a gauge field 
A 1 *, of the type C D can generate a bispectrum enhanced in the squeezed limit, 

with a specific signature in the form of a non-trivial dependence on the angle between the small 
and large wavevectors [121-123]. 

Large field models with an approximate shift symmetry also naturally contain a coupling between 
a pseudoscalar axion inflaton and a gauge field L D — (a/4/)</>.F /i! ,.F A “ / , where a is a dimensionless 
parameter and / is the axion decay constant. In these scenarios, gauge field quanta are produced 
by the background evolution of the inflaton, and in turn feed the curvature perturbation through 
an inverse decay process of the gauge field, resulting in a bispectrum shape strongly correlated 
with the equilateral template [124-128]. 

Isocurvature non-Gaussianities. — Multi-field inflationary models can generate residual pri¬ 
mordial isocurvature perturbations after inflation. Although their impacts on primordial power 
spectra are tightly constrained [14], they can contribute significantly at the level of the bispec¬ 
trum, producing in general both a pure isocurvature bispectrum and mixed bispectra because of 
their correlation with the adiabatic perturbation [129, 130]. 

8 Observational probes and constraints 

Non-Gaussianities can be constrained by studying the signatures they leave in the anisotropies 
of the CMB and in the large scale mass distribution in the Universe. Amongst the two, the CMB 
provides our cleanest window into the physics of the primordial universe, because its anisotropies 
are to a very good approximation linearly related to the primordial fluctuations. On the contrary, 
structure formation is an intrinsically non-linear process, and thus offers a more convoluted probe 
of the initial conditions. The CMB data, and in particular recent ones from the Planck mission, 
provide us with the best constraints so far [131, 132], and we briefly explain in the following how 
they have been derived (see the contribution by Bouchet to this volume for more details). We 
then turn more qualitatively to LSS, from where drastic improvements should originate in the 
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future [133]. 


Primordial non-Gaussianities in the CMB and Planck constraints. — The CMB maps 
of temperature and polarisation anisotropies are analysed by means of a decomposition into 
spherical harmonics: 

-fin) = ^2 a Jm Y em(n) 

lm 

E ip) = 

lm 

where the link between the CMB multipoles af m (X = T, E) and the primordial curvature 
perturbation ( is well known at linear order (see the contribution by Durrer to this volume). Using 
statistical rotational invariance, the CMB angular bispectrum, i.e. the three-point correlator of 
the aim s, can be written as 

( n X 1 n X 2 f . X 3 \ _ uX 1X2X3 (o o\ 

\ a hm 1 a l2m 2 a hirL 3 l ~ y m 1 m 2 m 3 u £ 1 e 2 l 3 V 5 - 0 ! 

where Qm[mlm 3 i s a known geometrical factor related to Wigner 3j-symbols, and bf^g X3 is re¬ 
ferred to as the reduced bispectrum. It is the counterpart in i space of the primordial shape S in 
Eq. (3.6), to which it is linearly related. Importantly, the signal one is looking for is way to small 
for a mode-by-nrode measurement of the (reduced) bispectrum. Instead, one compares theoret¬ 
ically motivated bispectrum templates to the observed one, and fit for their overall amplitudes 
(/nl’/nl The corresponding optimal estimator is well known [48], but its implementation 
is very challenging. The computational cost of its direct evaluation is totally prohibitive for 
high-resolution data like the ones of Planck, and several approaches have been developed to 
circumvent this problem. A central idea is to use templates that are separable, i.e. written as 
linear combinations of products of functions of respectively ki,k 2 ,k%, which renders the problem 
tractable [134]. The local shape (6.6) is of this type, while the equilateral and orthogonal ones 
(Eqs. (5.11)-(5.12)) are separable approximations of original non-separable shapes, and were de¬ 
vised with this data-analysis perspective in mind. Building on this idea, modal estimators are 
based on decomposing the bispectrum (both the theoretical ones and the observed one) into 
a sum of uncorrelated separable templates forming a complete basis in bispectrum space, and 
measuring the amplitude of each [135, 136]. This way, one can reconstruct once and for all the 
total bispectrum (see Fig. 5) and compare it to any theoretical shape, separable or not. At 
last, by exploiting the fact that angular bispectra of interest are generally smooth functions in 
harmonic space, binned estimators, which do not rely on separability, have also been developed 
and provide a complementary approach [137]. 

Based on these methods, the Planck collaboration performed a comprehensive analysis of 
their data, multiplying tests and cross-validations, quoting as their result for the three main 
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Figure 5. Planck modal reconstruction of the CMB temperature reduced bispectrum, plotted using 
several density contours [132]. 

shapes: 

( /nl = 2-5 ± 5.7 \ 

/n! = -16 ± 70 (68 % CL, T) 

VNL h = ~ 34 ± 33 / 

from temperature data alone, and 

( f%£ = 0-8 ± 5.0 \ 

= -4 ± 43 (68 % CL, T + 

V/NL h = " 26 ± 21 / 

from combined temperature and polarisation data. Despite having passed several consistency 
and robustness tests, it should be said that the polarisation data are known to suffer from some 
systematics and the (T + E) constraints should hence be considered as preliminary. 

Besides these estimates, the Planck collaboration constrained an impressively wide range 
of NG bispectra, including all the ones mentioned in this review and more. While they found 
no evidence for primordial non-Gaussianities, the ‘hints’ of NG reported in the Planck 2013 
analysis of oscillatory patterns are confirmed and reinforced with a more rigorous statistical 
analysis. To be precise, there appears to be no evidence for any individual feature or resonance 
model at a particular frequency, but multiple feature or resonance models might explain the 
apparently high non-Gaussian signal observed in the bispectrum reconstruction (see Fig. 5). The 
evidence of several of these models increased to more than 3cr with the additional polarisation 
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information, taking into account the ‘look-elsewhere’ effect, i.e. the fact that when scanning 
across the parameter space of models, Gaussian noise can lead by chance to a large apparent 
signal. No strong claim can be made at this stage though, as any information gathered from the 
polarisation data can only be considered as preliminary. 

At the level of the trispectrum, besides the constraint tnl < 2800 (95 % CL) from the Planck 
2013 analysis [131], the Planck 2014 analysis further gives ^nl = (—9.0 ± 7.7) x 10 4 , as well as 
estimates of two intrinsically quantum trispectra. 

Although the study of the angular bispectrum and trispectrum provides optimal constraints, 
let us point out that very good and consistent constraints on and ^nl were also obtained by 
using completely different tools, known as Minkowski Functionals, and which are based on the 
study of morphological features of random fields [138]. 

Eventually, one word about non-primordial NG (see the contribution by Bernardeau to this 
volume for further details): the high-quality data of Planck are such that, despite the small 
amplitude of primordial fluctuations, nonlinear effects from General Relativity itself, which are 
present even for perfectly Gaussian initial conditions, have to be taken into account. Most no¬ 
tably, the correlation between the gravitational lensing of the CMB anisotropies and the so-called 
integrated Sachs-Wolfe effect gives rise to a secondary CMB bispectrum that bias the naive esti¬ 
mate of /jJl by about 5. This contribution, of the same order as the statistical error in (8.4)-(8.5), 
has been subtracted from the results. 

Primordial non-Gaussianities in LSS. — An obvious way to probe primordial non-Gaussianities 
in LSS is via the study of higher-order correlators of the density field of galaxies, like the halo bis¬ 
pectrum. However, extracting the primordial information from such measurements is hampered 
by the significant non-linearities caused by the gravitational clustering itself. Hence, despite 
recent progress, we do not currently understand the galaxy density contrast to the accuracy 
required by the galaxy surveys, and there is no constraint to date on primordial NG from the 
halo bispectrum. 

There have been impressive progress in another direction though, called the scale-dependent 
bias. LSS surveys do not directly observe the mass distribution in the universe, but rather 
luminous biased tracers of it like the distribution of galaxies. In the presence of local type NG, 
it has been discovered recently that the bias between the density contrast of tracers and the one 
of matter acquires a non-trivial scale-dependence behaving like 1 /k 2 at large scales, making it 
a very sensitive probe [139]. Intuitively, NG in the squeezed limit implies that long-wavelength 
fluctuations in the gravitational potential (equivalently the curvature perturbation) locally rescale 
the amplitude of small scale matter fluctuations, thus affecting the threshold needed to form 
(halos of) galaxies, and thus the bias, the 1 /k 2 effect coming from the Poisson equation relating 
the Laplacian of the gravitational potential to the matter over-density. This effect manifests itself 
at the level of the halo power spectrum , which greatly simplifies the analysis compared to higher- 
order correlators. Using photometric quasars data from SDSS, the most stringent constraints to 
date give [140] —49 < < 31 and —2.7 x 10 5 < ^nl < 1.9 x 10 5 when constrained individually, 

and —105 < /j^ < 72 and —4.0 x 10° < <7nl < 4.9 x 10° in a joint analysis (all at 95 % CL). These 
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constraints are weaker than recent ones from the CMB, but provide a useful complementary and 
consistent picture. 

9 Implications for early Universe physics and perspectives 

The simplest single-field slow-roll models of inflation recently passed very stringent tests due 
to the lack of measurable deviations from Gaussianity of the primordial fluctuations. The con¬ 
straints that have been obtained strongly limit the different alternative mechanisms that have 
been proposed to explain the seeds of cosmological perturbations (see Ref. [132] for detailed 
constraints on numerous specific scenarios). 

The bounds on equilateral and orthogonal NG translate into a limit on the speed of sound 
of general single-field models of inflation (see section 5): 

c s > 0.020 95 % CL (T—only) (9.1) 

and c s > 0.024 (95 % CL) with the addition of the preliminary polarisation data. In the language 
of the effective field theory of inflation [141, 142], where inflation is studied as the theory of 
fluctuations of spontaneously broken time translations around a quasi de-Sitter background, the 
NG parameters can be related to the energy scale of the inflaton self-interactions as C/nl” ~ 
H 2 / A 2 (see [143-145] for details). Using these arguments, one can show that the cosmological 
dynamics is unaffected by higher-derivative terms for f^ OIth 1$ 1, corresponding to c s ~ 1. 
While current constraints still allow inflation to have wildly different dynamics than slow-roll, 
reaching the sensitivity A/^j]° rtl1 — 1 hence represents a well defined theoretical target for future 
experiments. 

From the single-field consistency relation (5.13), to local type non-Gaussianities in multifield 
inflation, through quasi-single-field inflation, we have seen that the behaviour of higher-order 
correlation functions in squeezed limits is a sensitive probe of the number of degrees of freedom 
active during inflation. In this respect, the stringent bound (8.4)-(8.5) obtained on severely 
constrains the dynamics of models in which light degrees of freedom beyond the inflaton play a 
role on super-Hubble scales. Compared to the bound on it is however harder to interpret 

it without specific mechanisms in mind, as the link between multifield scenarios and observables 
is very much model-dependent. In the large class of spectator models though, in which inflation 
is driven by a field while another one generates the primordial fluctuations, I/^lI ^ 0(1) is a 
rather generic prediction [146], like in the curvaton scenario or in modulated reheating. While 
achieving the sensitivity to test the single-fied consistency relation is currently out of reach, 
reaching A~ 1 therefore provides a challenging but conceivable target of particular interest. 

Due to the limitations set by Silk damping and foregrounds, it is unlikely that the CMB 
will offer significant improvements in the characterization of the statistical properties of the pri¬ 
mordial fluctuations. The impressive constrains on primordial non-Gaussianities from the Planck 
collaboration have indeed essentially achieved optimality, as shown by comparison with Fisher 
matrix forecasts. The way forward here seems to be LSS, including 21cm line radio surveys [147], 
which probes complementary smaller scales and enables one to extract genuine 3D information, 
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in contrast with the two-dimensional CMB sphere. Realistic optimised galaxy surveys have the 
potential to reach A~ 1 by using the scale-dependent galaxy bias [148]. As for other types 
of NG, one will have to resort to higher-order correlators like the halo bispectrum. Although one 
is hampered by theoretical uncertainties and observational systematics, the theoretical target 
A/^L° rth — 1 is within reach of currently planned surveys [133]. “When you have exhausted all 
possibilities, remember this: you haven’t.” supposedly said Thomas Edison. 
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